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We summarize a new analysis of the distribution <fc* of charged leptons produced in de- 
cays of Z and 7* bosons in the Collins-Soper-Sterman (CSS) formalism for transverse 
momentum resummation. By comparing the <j>* distribution measured at the Tevatron 
with the resummed CSS cross section with approximate 0(a^) Wilson coefficients, we 
constrain the magnitude of the nonperturbative Gaussian smearing factor and analyze 
its uncertainty caused by variations in scale parameters. We find excellent agreement be- 
tween the <f>* data and our theoretical prediction, provided by the ResBos resummation 
program. The nonperturbative factor that we obtained can be used to update resummed 
QCD predictions for precision measurements in inclusive W and Z production and for 
comparisons to various models of nonperturbative dynamics. 
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Leptonic angular distribution <j>* as a probe of nonperturbative dy- 
namics. Precision of the current data from the Tevatron and LHC on the inclusive 
transverse momentum (Qt) distributions of W and Z bosons imposes growing de- 
mands on theoretical calculations. Qt distributions are used for lucent tests of 
QCD factorization and in measurements of the W boson mass that place impor- 
tant constraints on the parameters of electroweak symmetry breaking. A variety of 
radiative contributions affect the Qt distribution of a heavy boson at the current 
level of accuracy. They include NNLO QCD and NLO EW perturbative corrections, 
logarithmic QCD contibutions that dominate da / cLQt when Qt — ► 0, and also non- 
perturbative power-suppressed terms that modify the Qt distribution when Qt is 
below a few GeV. Collins, Soper, and Sterman (CSS) [TJ [2j [3] have developed a 
QCD factorization approach to include all such terms order by order in a s . The 
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CSS formalism combines the resummation of large Sudakov factors in the small-Qy 
limit [H O |6] with fixed-order QCD contributions at large Qt- As a recent devel- 
opment, the NNLL/NNLO expression for the resummed Qt distribution has been 
published [7J. The CSS method is a classical realization of the factorization based 
on transverse-momentum dependent (TMD) distribution functions, in which PDFs 
and fragmentation functions depend explicitly on intrinsic transverse momentum 
in addition to the usual momentum fraction variables. While theoretical methods 
of TMD factorization [H El [113 HH [HI H31 HH Ull [13 HH US] and soft-collinear 
effective theory HOI IS1 (HI [2H H2] undergo rapid developments, the CSS formalism 
is well-suited for detailed phenomenological studies, as it is implemented in detail in 
practical simulations. In this short paper, we summarize a recent application of the 
CSS formalism for the computation of the angular distribution da/d<f>* of the lepton 
pairs in Z/^y* boson production at the Tevatron. Full details will be presented in 
a separate forthcoming paper [24]. The <j>* distribution is closely related to the Qt 
distribution in the limit </>* oc Qt /Q 0. Its analysis may provide valuable insights 
about the nonperturbative QCD dynamics. 

The D0 collaboration has published detailed measurements of the cj>* depen- 
dence in the electron and muon decay channels [25]. We ask if the D0 data cor- 
roborate the universal behavior of the resummed nonperturbative terms that was 
observed in the global analyses of Qt distributions in 7* and Z production at fixed- 
target and collider energies [13 [57] . We also investigate the rapidity dependence of 
the nonperturbative terms, which may be indicative of new types of higher-order 
contributions |28j . The small-QT part of the distribution is obtained in the CSS for- 
malism by the Fourier-Bessel transform of a form factor W(b, Q) that depends on 
the transverse position variable b. For large boson virtualities Q of order 100 GeV, 
the overall form of the Qt distribution at all Qt is determined by small-fe contri- 
butions, arising at energy scales /j~ l/!>> 1 GeV. Such contributions are entirely 
predicted in perturbative QCD theory. They dominate the cross section at any Qt 
value [5[ [29] . When Qt is below 5 GeV, the production rate is also mildly sensitive 
to the behavior of W(b, Q) at b > 0.5 GeV -1 , where the perturbative expansion is 
increasingly unreliable because of the vicinity of the Landau pole in a s {l/b). We 
are interested to know which forms of the large-6 extrapolation of W(b, Q) are com- 
patible with the observed behavior of Qt and <fi* distributions. A comprehensive 
solution for W(b, Q) at large b remains elusive, but efforts to derive it produced 
several instructive models, such as [30] [31] [32] [33J [34l [35J [20] . From the phenomeno- 
logical point of view, a typical Z data set does not have sensitivity to distinguish 
between these models. In deeply nonperturbative region that is characterized by b > 
1 GeV -1 , the form factor W(b,Q) is strongly suppressed and does not contribute 
to da/dQT for Q of order Mz [27] . Only the contributions from the transition 
region of b of about 1 GeV -1 , where fiKi 1/6 is about 1 GeV, are numerically non- 
negligible compared to the leading-power cross section predicted by perturbative 
QCD theory. In the transition region, the extrapolation of the perturbative expres- 
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sion provides a reasonable approximation for the leading-power (logarithmic) part 
of W(b,Q). It can be realized in the "revised fr* model" [27], which parametrizes 
the extrapolated contribution by a flexible form that depends on a single param- 
eter b max . In addition, a few suppressed terms proportional to even powers of b 
play some role in this interval of b. In comparisons to the experimental data, their 
cumulative effect can be usually approximated by a single Gaussian smearing fac- 
tor exp{— a(Q)b 2 }, where the coefficient a(Q) is found from the experiment. This 
arrangement provides a few-parameter approximation for viable nonperturbative 
models in the phenomenologically relevant region of b. In our previous work [27) . 
the magnitude and Q-dependence of the Gaussian factor were determined from Qt 
distributions of the Drell-Yan pairs. Recently, the 0* distribution has been proposed 
as a sensitive probe of the small-Qr dynamics [36] . as it has reduced uncertainties 
associated with the lepton momentum resolution. Here we update the constraints 
on the Gaussian smearing factor in Z boson production using the 4>* distribution. 

Our analysis is carried out using the program ResBos [37l ES [38] , which real- 
izes the CSS formalism to compute fully differential cross sections of lepton pairs 
in production of high- mass virtual photons (7*) and heavy electroweak bosons (W, 
Z, and H). A resummed treatment of new variables <It and 4>* and their relation- 
ship to Qt was studied in (39] 00]. The fully differential output from ResBos can 
also be cast in the form of the (f>* distribution. In the new analysis, we find an 
excellent agreement between the ResBos prediction and the cj>* data, contrary to 
the conclusion made by the D0 paper [25) . However, the quality of agreement de- 
pends on the inclusion of perturbative loop contributions and selection of scales in 
the resummed cross section. In Refs. [3^1 SOI SU HH H]> the evidence for a nonzero 
nonperturbative factor in Z production has been contested in the light of the uncer- 
tainty in the resummed cross section due to the dependence on factorization scales. 
The evidence for nonperturbative smearing is inconclusive at the NLL+NLO level 
because of a large scale dependence. To address this point, we fully include the 
scale dependence in the small-Qr part of the resummed cross section up to 0{a 2 s ) 
[24] . Despite the scale uncertainties, the statistical analysis of the fit to the (p* data 
indicates pronounced preference for a nonperturbative Gaussian contribution with 
a(Mz) PS 1 GeV -2 . The main reason is that the Gaussian suppression of the large-6 
tail of W(b,Q) alters the resummed da/dQr in a different way than the factor- 
ization scales in the leading-power part of W(b, Q). The nonperturbative Gaussian 
factor suppresses the rate only at Qt below 2-3 GeV, while the leading-power scale 
dependence affects a broader interval of Qt values. This results in a characteristic 
shift of the peak in the da /cLQt distribution due to the nonperturbative suppression, 
which is distinct from the typical scale dependence. 

Particulars of the calculation. In the CSS formalism, the full resummed 
(RES) Qt distribution is commonly represented as a combination of the resummed 
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(W), finite-order (FO), and asymptotic (ASY) terms [3]: 

do \ ( do \ I da \ /do 



dQr J res \dQrp J yy \dQrp J pQ X^-Qt / asy 



(1) 



The accuracy of the DO measurement demands that the main perturbative correc- 
tions due to the QCD and electroweak radiation are included. Numerical, but not 
analytical, results for the complete NNLO QCD contribution to the resummed Qt 
distribution in Z/j* production at the Tevatron have been recently released [7]. 
ResBos includes the dominant NNLO QCD contributions and provides a faithful 
estimate for the remaining small NNLO contribution that has not been published in 
an analytical form, as summarized below. Thus, effectively ResBos is close to the 
full NNLO precision in the kinematical region relevant for this analysis. The elec- 
troweak (EW) corrections to Z production compete in magnitude with NNLO QCD 
contributions and are available to NLO [43j U4j |45l 06] . In the comparison to the DO 
data, we do not include the EW corrections in our theory prediction, but correct 
the fitted 0* data by subtracting the predominant correction due to the final-state 
photon radiation obtained by the Photos code [47]. Upon the inclusion of these 
contributions, scale dependence remains a major systematic uncertainty affecting 
our theory prediction. In the latest ResBos implementation, the fully differential 
CSS resummed cross section is given by [3 7) 

do (/ii/i 2 -> (Z ->• t£)X) _ 1 Q 2 

dQ 2 dydQ 2 T dVl ~~ 48irS (Q 2 - M 2 ) 2 + Q 4 T 2 Z /M 2 

0e« ]T Wf ert (K,Q,x 1 , X2 ,Q,C 1 ,C 2 ,C 3 )W NP (b,Q) 

j—u,d,s... 

+Y(Q T ,Q,x 1 ,x 2 ,il,C^, (2) 

in notations of Ref. [37 . The leading-power ( "perturbative" ) form factor W pert is 
defined by 



WPert = \H(Q,n,C±)\ 2 exp 



f 2Q ^AfcCJlnl^Pj +B( } l;C 1 ,C 2 ) 



x £ [c 3a ® f a/hl ] L, <f \ ® h/h*] U, % y ) (: *' 

a=g,q \ 2 / b=g q 



in terms of the hard part H(Q, fi, C4) (dependent on Q and the solid angle ft = 
of Z boson decay in the Collins-Soper frame), Sudakov exponent, and con- 
volutions \Cj/ a ® fa/n\ 01 Wilson coefficient functions Cj/ a (z, Ci/C 2 , Hf — C3/6) 
and parton distribution function f a /H{ z i A*f) for a parton a inside the initial-state 
hadron h. W NP (b,Q) is the nonperturbative factor explained below. Y is the dif- 
ference between the finite-order and asymptotic cross sections and it dominates 
at Qt of order Q. Constants Cj (for i = 1,...,4) are the scale coefficients that 
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determine several factorization scales introduced by resummation. They arise be- 
cause W pert depends on two distinct momentum scales 1/6 and Q. Specifically, 
C\ = bfi and C 2 = n/Q, where fj, is the scale introduced by the evolution equa- 
tions that control the large logarithms. The constant C3 arises in the Wilson 
covolutions[Cj / ® f a /H\ an d specifies the factorization scale (j.f = C 3 /b in the 
PDFs f a /H{ z T Hf)- C4 = hh/Q is the scale constant in the hard part H of the 
resummed term W and the Y piece, which we select as C4 = C 2 for simplicity. 
The expression for the resummed cross section becomes particularly simple for a 
"canonical" combination of the scale coefficients, given by C\ = bo, C 2 = 1, C3 = bo, 
where bo — 2e~~ tE , and je = 0.577... is the Euler-Mascheroni constant. We eval- 
uate Y to 0(a 2 ) based on the calculation in [48j |49l [50] . The functions A(p; C\) 
and B{ji)Cx,C%) are computed up to orders 0(a^) and 0(a 2 ), respectively. The 
Wilson coefficient functions are computed exactly up to 0(a s ). The only unavail- 
able part of the NNLO resummed cross section is the 0(a 2 s ) Wilson coefficient C^ 2 } 
which receives contributions from two loop virtual diagrams. However, from the 
fixed order NNLO calculation [50] this contribution is small in magnitude (2-3%), 
mostly affects the overall normalization of the W term, and has weak kinematical 
dependence. Thus, without losing accuracy, one can approximate this term by 

cf}{z, Cx/C 2 ,C 3 ) = 6C& 6(1 - z) 6 ja + L(C U C 2 , C 3 ), (4) 

where L(C\,C2tC^) = for the canonical combination. Since 5C^ is nearly con- 
stant in the kinematical region relevant for Z production, we estimate its magnitude 
from the known value of the 0(a 2 s ) if- factor for the inclusive cross section da/dQ 
that is known for a long time [51] and was evaluated in our analysis by the computer 
code Candia [52J [53] . The inclusion of the estimated SC^ in the calculation has 
practically no effect on our conclusions. Finally, L(Ci, 62,63) is found exactly by 
requiring the independence of the a s series expansion of W on the choice of C\, 
C2, and C3 order by order. By truncating the series at 0(a 2 ), we must have the 
equality W(b,Q,C 1 ,C 2 ,C 3 )\ 0{a 2 ) = W(b,Q,Ci = C 3 = b ,C 2 = l)|o( a J), which 
allows us to completely reconstruct the dependence of the 0(a 2 ) part on the scale 
parameters Ci. The dependence on the scale constants is illustrated in Fig.Q] which 
shows the ratio of the experimental and best-fit theoretical values of (l/<x) ■ da/d<f>* 
for electrons with a constraint \yz\ < 1 on the Z boson rapidity. The best agreement 
with the data are obtained for {Ci = 2&o, C 2 = 1/2, C3 = 26o}- Several error bands 
are obtained by variations of the indicated scale parameters around this best-fit 
combination. The variation of the scales affects the quality of the fit quantified by 
X 2 and modifies the cross section in a wide range of <fi*. In contrast, the variation 
of the nonperturbative factor is pronounced only at 0* < 0.5, which corresponds 
to typical Qt of a few GeV. This difference allows one to discriminate between the 
nonperturbative Qt smearing and perturbative scale dependence. 

The nonperturbative factor. In order to extrapolate W pert in Eq. ([2]) to the 
large b values of order or above 1 GeV -1 , we evaluate it in the revised 6* (b, b rnax ) 
model [27] as a function of 6* = b/y/T+ {b/b max ) 2 with b max — 1.5 GeV -1 . 
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Fig. 1 . Impact of the scale variation on the agreement between theory and data for \yz\ < 1 ■ 

The factorization scale in the convolutions [C (8> /] is set in this model to fip — 
C l 3/6*(&, bo/Qo), where Qo = 1 GeV is the initial scale of the PDFs. This choice of 
bmax is preferred by the global fit to Drell-Yan Qt data, where it both improves the 
agreement with the data and preserves the exact form of the perturbative expansion 
for W pert at b < 1 GeV -1 . The (f>* data is dominated by the narrow vicinity of Q 
around Mz- Therefore, for the purpose of this analysis, it sufhcies to approximate 
the nonperturbative Gaussian factor as 

W NP (b, Q a M z ) = cxp [-Oi(Afe) b 2 ] . (5) 

This expression simplifies a more general parametrization [271 126) 



W JslJ {b,Q) = exp 



(6) 



in which the coefficients a±, 02, and 03 can be separated by fitting to several data 
sets at distinct y/s and Q combinations. Once a\(Mz) is known, the coefficients 02 
and 03 from the global Qt fit of Ref. [27] can be used to find W NP in W boson 
production. In each yz bin of the electron and muon </)* data, we compute \ 2 an d 
use it to determine the nonperturbative coefficient a\. The results are shown in 
Fig. [5J In order to estimate the impact of the scale dependence, the fit includes a 
correlation matrix quantifying the uncertainty in the theory cross section due to 
variations in the scale constants C\,Ci and C3. In this case, we assign the 68% 
confidence to the interval — 1 < \og 2 (Ci/C^ est ~^ t ) < 1. According to the figure, all 
rapidity bins generally prefer a non-zero a±, even if the scale shifts are included. 
The scale dependence increases the errors and makes them very asymmetric, but 
a downward change in the best-fit a\ is more disfavored than the upward change. 
With the scale shifts, we obtain the central value of a±(Mz) = I.OS^q'^ GeV" 2 . No 



September 11, 2012 1:52 WSPC/INSTRUCTION FILE qcde2012v4 



Z/7* Transverse momentum distribution at hadron colliders 7 



> 
a 
O 





Without scale shifts 








■ e j 







l.vzl<l 



l<l)'zl<2 
yz bins 



l.vzl>2 




\y 2 \<\ K|y z |<2 \y z \>2 

y z bins 



Fig. 2. Estimate of coefficient ai in the three bins of rapidity of the data without (left inset) and 
with (right inset) C\,Ci and C3 shifts. 



significant dependence of the best-fit a\{Mz) values on the rapidity yz is observed, 
but the uncertainty is very large in the \yz\ > 2 bin. 

Conclusions. A new version of the resummation program ResBos was em- 
ployed to examine the differential cross section at small values of the leptonic angle 
<fi* at the Tevatron. The data on the Z/-f* (f>* distribution collected by the D0 
collaboraton was compared with the CSS resummed cross sections with approx- 
imate 0{a 2 s ) Wilson coefficient functions and complete 0{a 2 s ) scale dependence. 
ResBos agrees well with these data, provided we use the hard scales Q/2 (i.e., the 
scale coefficients C2 = C4 = 1/2) in the resummed cross section. We determined 
the nonperturbative factor W NP (b, Q) preferred by the 0* , while evaluating the ef- 
fects that have comparable magnitude: QCD scale dependence at 0(a 2 ) and NLO 
electromagnetic contributions. In a fit that allowed for variations of the resumma- 
tion scale parameters Ci, we observe a 2.5a preference for a non-zero Gaussian 
smearing factor a\{Mz) ~ 1.1 GeV~ 2 . This central value is in agreement with 
the findings of the previous analyses of Z Q^-dependent distributions [221 133 HZ] • 
The non-pertubative factor W NP (b, Q) is a part of the complete resummed factor 
W pert (b*, Q)W NP (b, Q), which is well-constrained in the phenomenologically rele- 
vant region b < 1 GeV -1 by a combination of the 0(a 2 s ) PQCD calculation and the 
nonperturbative factor found from the experimental data. Further information on 
this study will be provided in Ref. [24! . It will be used to update ResBos predictions 
for future Tevatron and LHC studies. 
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